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Context. The cosmic time around the z~l redshift range appears crucial in the cluster and galaxy evolution, since it is 
probably the epoch of the first mature galaxy clusters. Our knowledge of the properties of the galaxy populations in 
these clusters is limited because only a handful of z~l clusters are presently known. 

Aims. In this framework, we report the discovery of a z~0.87 cluster and study its properties at various wavelengths. 
Methods. We gathered X-ray and optical data (imaging and spectroscopy), and near and far infrared data (imaging) in 
order to confirm the cluster nature of our candidate, to determine its dynamical state, and to give insight on its galaxy 
population evolution. 

Results. Our candidate structure appears to be a massive z~0.87 dynamically young cluster with an atypically high 
X-ray temperature as compared to its X-ray luminosity. It exhibits a significant percentage (~ 90%) of galaxies that 
| are also detected in the 24/xm band. 

Conclusions. The cluster RXJ1257. 2+4738 appears to be still in the process of collapsing. Its relatively high temperature 
is probably the consequence of significant energy input into the intracluster medium besides the regular gravitational 
infall contribution. A significant part of its galaxies are red objects that are probably dusty with on-going star formation. 
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1. Introduction 
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On the observational side, we might expect selection ef- 
fects to favor finding the most X-ray luminous clusters when 
searching for clusters at high rcdshift. This is not the case, 
however. For example, the compilation by Ota et al. (2006) 
shows that 50% of the cluster sample between z = 0.3 and 
0.56 have a bolometric X-ray luminosity (Lx,boi) greater 
than 10 45 ergs s _1 compared to only 20% of the sample 
above z = 1 (Fassbender 2008, Kocevski et al. 2009). These 
X-ray luminosities suggest that beyond z=l, clusters are 
growing, but have not reached their peak X-ray luminosity. 

In comparison, we estimate the free fall times for clus- 
ters (with a typical radius of 1 Mpc and masses between 
(0.5 - 5) x 10 14 Mq) of the order of 1-2 Gyr (Sarazin 
1986). The elapsed time between z = 2 and z = 1 is about 
2.6 Gyrs. Thus, at redshift range near ~ 1 we can expect 
mainly to find clusters of (0.5 — 5) x 10 14 M© that have just 
completed their initial infall (see Sadeh & Rephaeli 2008 
and references therein). 

Since X-ray detectable clusters near z = 1 are probably 
just forming large enough systems to become detectable, a 
study of the relationship between the galaxy population and 
the X-ray emission is particularly instructive for refining 
models of both large scale structure evolution and galaxy 
evolution. See Ettori et al. 2004b, Andreon 2006, Ota et al. 
2006, for comprehensive studies of X-ray detected clusters 
and cluster evolution, and Cooper et al. 2008, Yan et al. 
2008 for recent studies on galaxy formation and evolution 
related to clusters of galaxies at z~ 0.8 — 0.9. 

In the process of our search for moderately distant (i.e., 
z ~ 1) clusters of galaxies that are detectable via their X- 
ray emission (Adami et al. 2000, Romer et al. 2000, Ulmer 
et al. 2005, Adami et al. 2007), we report here the discovery 
of a cluster with z = 0.866, RX J1257.2+4738 (hereafter 
for brevity referred to as RX J 1257). The cluster was first 
detected in X-rays in the ROSAT data archive, and the i'- 
band and Ks-band follow-up found a concentration of red 
galaxies that could be explained by the existence of a cluster 
with z <> 0.6. 

This cluster is particularly interesting because it con- 
tains several red galaxies with [Oil] 3727 emission as well 
as several others that are detected both in the Spitzer MIPS 
24 /jm and Spitzer IRAC bands. Without the detected 
emission lines or the MIPS emission, these red galaxies 
might have been characterized as early type galaxies, or 
red-dead galaxies (Elbaz et al. 2007). The emission lines im- 
ply young galaxies undergoing star formation and the high 
MIPS flux implies dusty late type galaxies with ongoing 
star formation. We report here a detailed study based on 
extensive observations including Chandra, XMM-Newton, 
Spitzer, Gemini, Subaru, and ARC. Also, as a by-product, 
we found candidates for lensed (and hence magnified) back- 
ground galaxies or AGNs that are at least at z > 3, and 
could be at redshifts as high as 10. These will be presented 
in a separate later paper. A preliminary version of the re- 
sults reported here were given by Ulmer et al (2008). 

Throughout, we assume a concordant ACDM cosmology 
with H =71 km s" 1 Mpc" 1 , fi A ,o = 0.73, ft m ,o = 0.27, from 
which we compute the angular scale as 7.72 kpc/arcsec, a 
luminosity distance (Dl) = 5544.7 Mpc and an angular 
distance (D A ) = 1592.4 Mpc (Wright 2006). 



2. Observations, Analysis and Results 

2.1. X-rays 

The net Chandra exposure time was 37 ksec in very faint 
mode with ACIS-S3 (back illuminated chip) on July 24, 
2007. 

The XMM-Newton net exposure time, after high back- 
ground filtering was 16 ksec (11 and 5 ksec for each usable 
exposure). The nominal exposure time (83 ksec allocated) 
was heavily polluted by flares. The observations were made 
with the Thin 1 filter in Prime full window on May 24th 
through 26th, 2007. 

2.1.1. Chandra data processing 

The data were reduced with the CIAO version 3.0.10 fol- 
lowing the Standard Data Processing, producing new level 
1 and 2 event files. 

We further filtered the level 2 event file, keeping only 
ASCA grades 0, 2, 3, 4 and 6, and restricted our data re- 
duction and analysis to the back-illuminated chip, ACIS-S3. 
We checked that no afterglow was present and applied the 
Good Time Intervals (GTI) supplied by the pipeline. We 
then checked for flares in the light-curve in the [10-12 keV] 
band. No flares were detected. 

2.1.2. XMM-Newton data processing 

The XMM-Newton data were so badly contaminated with 
flare events, that in addition to the standard processing 
tools, visual inspection of the overall count rate was needed 
as well in order to select only good data. 

RX J1257 was observed in standard Full Frame mode 
with the "thin" filter with the two EPIC MOS1 and MOS2 
and the PN detectors. The basic data processing (the 
"pipeline" removal of bad pixels, electronic noise and cor- 
rection for charge transfer losses) was done with package 
SAS V5.3, thus creating calibrated event files for each de- 
tector. 

For the MOS1 and MOS2 cameras, following the stan- 
dard procedure, we have discarded the events with FLAG ^ 
and PATTERN > 12; for the PN, we have restricted the 
analysis to the events with PATTERN < 4 and Flag =0. 

With the cleaned event files, we have created the re- 
distribution matrix file (RMF) and ancillary response file 
(ARF) with the SAS tasks rmfgen and arfgen for each 
camera and for each region that we have analyzed. 

2.1.3. Imaging 

We detect the diffuse emission of the cluster both in 
Chandra and XMM-Newton observations. If we model the 
XMM-Newton emission with an ellipse it has a semi major 
axis of 1.2 arcmin and a semi-minor axis of 1.0 arcmin (556 
kpc and 463 kpc at z = 0.866). The major axis runs along 
the E-W direction and the cluster is centered at 12h 57m 
12.2s, +47° 38' 06.5" (J2000). A point source we detected 
with Chandra appears embedded inside the diffuse emis- 
sion. We will discuss this source in § 12.1.41 and in § 13.2.11 
below. 



http : //asc .harvard. edu/ ciao/ 
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The X-ray emission is statically significant. The count 
rate of the Chandra observation, inside the ellipse where 
the cluster spectrum is extracted (see below), is 0.0065 ± 
0.0002 cnt/s/arcmin 2 (where the point sources are already 
masked out). On the same CCD, we estimate the back- 
ground count rate as 0.00495 ± 0.00015 cnt/s/arcmin 2 . 
Therefore, the net count rate of the cluster in the central 
1.2' x 1' region is 0.0016±0.0004 cnt/s/arcmin 2 , or 4er above 
zero. Combining all XMM observations (all detectors and 
exposures), we obtain for the same region (masking the 
point sources taking into account the larger XMM PSF) a 
net count rate of 0.0059±0.0002 cnt/s/arcmin 2 , well above 
a null count rate. 

We are dealing here with a densely populated field of 
active objects, which dominates the X-ray emission of the 
area. We therefore masked these point source contributions 
to determine the diffuse emission of the cluster. The ability 
to detect and verify the point-like nature of these objects 
was impossible to do using only XMM-Newton data due to 
a lack of angular resolution. We therefore used the Chandra 
data. The point source masks removed over 99% of the en- 
circled point source flux and thus the remaining flux was 
point source contaminated at a level well below statisti- 
cal uncertainty in our derived best fit spectral shapes and 
overall luminosity. 

For all the relevant masked X-ray point sources, we 
found optical counterparts that were bright enough to apply 
LePharqj. For those objects classified as galaxies, we esti- 
mated the probability that they belong to the cluster based 
on the probability distribution function of the photometric 
redshifts (hereafter the PDF). We determined the PDFs 
with the LePhare software The PDF gives the probability 
for an object to be at a given redshift. The probability for 
an object to be in a redshift interval is therefore obtained by 
integration over this interval and normalization to the total 
probability. We give below the list of the prominent X-ray 
point sources, the probability that they are in the cluster 
redshift interval, and the LePhare spectral classification. 

The main point sources showing up both in Chandra 
and XMM-Newton data are: 

- (a=194.2680, £=47.6377) associated with a bright galaxy. 
The probability of this bright X-ray object to be in the 
[0.82,0.90] redshift interval (potentially a cluster member) 
is lower than 0.2%. This is probably a late type active 
galaxy (type 19 from the LePhare spectral data base) 
equivalent to a starburst; 

- (a=194.2930, (5=47.6392) associated with a galaxy. The 
probability of this object to be in the [0.82,0.90] redshift 
interval is lower than 0.4%. This is probably a very active 
object, LePhare gives a better \ 2 with a QSO template 
than with a galaxy template; 

- (a=194.2930, (5=47.6339) associated with a blend of 3 
objects within a 2" long region. This object is discussed in 
detail in §3.3.; 

- (a=194.2570, (5=47.6251) associated with an active 
galaxy at z=0.662 from our Gemini spectroscopic catalog. 
This galaxy exhibits [Oil] and strong [OIII] lines, as well 
as broad H7 and H/3 emissions, characteristic of an active 
object; 

- (a=194.2950, (5=47.6476) which is probably a local star 
(not in the Simbad data base), since LePhare gives the best 

2 S. Arnouts & O. Ilbert; 
www.oamp.fr/people/arnouts/LE_PHARE.html 



X with a star template. Using the probability distribution 
function given by a galaxy template, the probability that 
the object's redshift is in the [0.82,0.90] redshift interval is 
lower than 9%; 

- (a=194.3051, (5=47.6132) associated with the GMOS 
guiding object (i' SDSS mag = 12.39). We observed its 
spectrum at the Observatoire de Haute Provence 1.93m 
telescope with the CARELEC spectrograph. The spectrum 
is typical of a local old star; 

- (a=194.2914, (5=47.6618) associated with a faint proba- 
bly high redshift object (from its PDF) which only has a 
probability of less than 1.5% to be part of the cluster; 

- (a=194.2876, (5=47.6641) associated with a faint proba- 
bly high redshift object (from its PDF) which only has a 
probability of less than 0.1% of being a cluster member; 

- (a=194.3070, (5=47.6715) associated with a relatively 
bright optical object. This object has a non-negligible prob- 
ability of being at the cluster redshift (12%) but since it is 
located at a projected distance of over 1 Mpc from the clus- 
ter center, it is probably not directly related with the struc- 
ture itself. We cannot exclude, however, that the galaxy is 
undergoing cluster infall. 

2.1.4. X-ray spectroscopy of the diffuse source 

The background spectra were extracted from the same ob- 
servations, on the same CCD chip, from a region close to 
the cluster, free of point sources. In this way, we took into 
account any count rate anomaly in the background. The 
center of the background region is only about 3.5 arcmin 
from the source, making the vignetting effect (for the cos- 
mic X-ray background term) negligible. The redistribution 
and ancillary files (RMF and ARF) were created with the 
SAS tasks rmfgen and arfgen for each exposure, taking into 
account the extended nature of the source. 

We have fit simultaneously spectra from the two PN 
exposures with the MEKAL spectrum model of a thermal 
plasma. The galactic neutral hydrogen along the line of 
sight was taken into account with the "phabs" photoelec- 
tric absorption model of Balucinska-Church & McCammon 
(1992). The hydrogen column density, Nh, was fixed to 
the galactic value given by the Leiden/Argentine/Bonn 
(LAB) Survey of Galactic HI (Kalberla et al. 2005); N H = 
1.2 x 10 20 cm' 2 . We have also fixed the metal abundance 
to Z = O.4Z and the redshift to z = 0.866 (see the optical 
spectroscopy section). In order to minimize any discrepancy 
between Chandra and XMM-Newton at low and high ener- 
gies, we limited the spectral fit to the 0.7-7.0 keV band. 

Since the X-ray field in the direction of the cluster is 
heavily contaminated by X-ray point sources, it is useful 
to archive the masking out of these sources as shown in 
Fig. [TJ We show the resulting best fit (3 model in Fig. [5] 
based on the assumption that the cluster has an intrinsi- 
cally smooth elliptical surface brightness distribution. We 
also show in Fig. [3] the 50% and 90% encircled energy radii 
to demonstrate that the X-ray emission is clearly diffuse. 
The assumption of a symmetric (3 model is only necessary 
for extracting the spectrum and was not used to determine 
the true shape of the X-ray emission. A detailed discus- 
sion of the nature of the inferences based on the shape of 
the emission is given below. We conclude, though, that we 
have have indeed detected diffuse emission. The results of 
the spectral fitting are given in Table [TJ 
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Fig. 1. X-ray imaging of CI 1257 in the 0.5-5.0 keV 
band, (a) Adaptivcly smoothed ACIS-S image with a mini- 
mum signal-to-noise ratio of 3 per smooth beam done with 
the COAI 3.4 csmooth task; (b) ACIS-S image smoothed 
with a fixed width Gaussian (kernel radius of 10 arcsec); 
(c) best-fit 2D /3-model; (d) composite image with all avail- 
able XMM data smoothed with a Gaussian kernel (radius 
6.4 arcsec). The green ellipses correspond to logarithmically 
spaced isocontours of the best-fit /3-model. The magenta 
circles show the masked regions (omitted both in the 2D- 
fit and in the spectral fits) corresponding to point-sources 
detected with Chandra data. The color bars indicate the 
imaging scaling in counts per pixel. 

Table 1. Summary of XMM-Newton Spectral Analysis 



kT 

(keV) 


f x (0.5-2.0 keV) 
(10 -14 erg s-'em -2 ) 


Lx,bol 

(hfolO^ergs- 1 ) 


xVdoff 


3.6 + _H 


2.2 ± 1.1 


2 +2 ' 9 


64.06/64 



' Best fit thermal bremsstrahlung model; Z = OAZq fixed. 



For future reference in the text, we use the spectral fit 
and the formula from Hilton et al. (2007) that is based on 
the temperature of the ICM, to derive R v =1.05 Mpc or 
2.3 arcmin, where R v is the virial radius of the cluster. 

Given the relatively low number of photons included 
in the spectrum (^1250 counts including background), the 
spectrum alone does not allow us to discriminate at a sta- 
tistically significant level between a thermal and power law 
spectrum that might be expected from an AGN. 

2.2. Visible, NIR and IR data 

In order to determine the nature of the galaxies in the clus- 
ter as well as to measure the redshift of the cluster, we 



10 100 
radius (physical pixels) 

Fig. 2. Radial surface brightness profile best fit /3-model 
to Chandra data after point source masking. The vertical 
dashed lines indicate the radii enclosing 50% and 90% of the 
energy from a point source, spread by the Chandra PSF. 



obtained visible, near IR and IR observations. We summa- 
rize the observations and the results of those observations 
below. 

2.2.1. Spitzer data reduction 

Spitzer (Werner et al. 2004) images obtained with the 
Infrared Array Camera (IRAC Fazio et al. 2004) were ob- 
tained on 2006 Dec 29. The observation consisted of 100- 
second integrations and 20 dithers giving a net exposure of 
2000 seconds in each of the four channels centered at 3.6, 
4.5, 5.8, and 8 /im and covering a ~ 5/2 x 5/2 field. A 
deep image at 24 /jm was made on 2007 June 02 with the 
Multiband Imaging Photometer for Spitzer (MIPS, Rieke et 
al. 2004) using 30-second integrations with a total exposure 
time of 9,050 seconds. The calibrated mosaic images from 
the standard Spitzer pipeline showed satisfactory removal 
of instrumental artifacts and were sufficiently clean for the 
IRAC data analysis. However, the MIPS 24 //m pipeline 
mosaic exhibited artifacts and gradients, which were sig- 
nificantly reduced through a zodiacal light subtraction, a 
self-calibrating flat, and an improved overlap correction al- 
gorithm. Source fluxes were obtained from these mosaic 
images with SExtractor (Bertin & Arnouts 1996) aperture 
photometry. 

The original units of the Spitzer images (MJy/sr) 
were translated to Jy for detected objects taking into ac- 
count the pixel size to which these images were rebinned 
(0.15x0.15 arcsec 2 ) and the general directions given in 
http://ssc.spitzer.caltech.edu/ □ 



Magnitudes were then computed following the same di- 
rections and translated from the Vega system to the AB 
system for the IRAC bands with the numbers listed in 
http://spider.ipac.caltech.edu/staff/gilli an/cal.html] 
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2.2.2. Gemini spectroscopy and imaging 

The measurement of redshifts in a z~0.9 galaxy structure 
required spectroscopic observations of galaxies as faint as 
i' ~ 23. In order to achieve this goal, we acquired 3 Gemini 
north GMOS observations (GN-2005A-Q-9, GN-2006A-Q- 
4, and GN-2006B-Q-38). We measured 45 redshifts of galax- 
ies with magnitudes between i' = 20 and 22.6 and exposure 
times varying between 3 and 4 hours. At the same time, we 
took advantage of the pre-image observing to acquire 2 deep 
optical images of this cluster in the i' and z' bands. The 
completeness limit of these images was close to i' —25 and 
z' —23.5 as shown in Fig. 0] (together with the completeness 
histograms in the J and Ks bands for comparison). 

We obtained images in the i' band under photometric 
conditions. In order to assess the quality of these magni- 
tudes, we compared our values with the SDSS estimates. 
SDSS images cover the same area with the same filters and 
are much less deep than our data. However, the brightest 
objects of our field are also detected in SDSS images. Based 
on 51 objects in common, we computed the mean value of 
the difference between SDSS i' band magnitudes and our 
estimates, and found a good agreement: 

i'sDSs ' ^'gmos = — 0.05± 0.17 

Images in the z' band were observed under non- 
photometric conditions. We therefore used the SDSS mag- 
nitudes to rescale our own measurements with 45 objects in 
common. The resulting dispersion between SDSS and our 
estimates is ±0.21 magnitude. 

We show in Fig. [3] the color magnitude relation based 
on the i' and z' data. Note that most if not all the scatter 
in the i'— z' direction is attributable to the dispersion in the 
derived values of z'. The cluster being at z— 0.87, the 4000 A 
break is included in the i' filter and located just before the 
z' filter. We also superpose on Fig.[3Jthe galaxies inside and 
outside the z=[0. 850, 0.874] range. Fig. [3] shows a possible 
(given the dispersion in z') concentration close to i'-z' —0.6 
(the red diamonds in Fig. [3]). 

We will argue in the following that our cluster is a dy- 
namically young structure, and thus that many of the galax- 
ies have probably not turned into ellipticals or SOs. The av- 
erage color of i' — z'=0.6 that we found is consistent with 
those of early spirals at z = 0.8 based on Fukugita et al. 
(1995) who predicted a color of —0.6 for these objects. 

Fig. [3] also exhibits a prominent concentration of bright 
and probably nearby galaxies with i'—z' —0.3. We retrieved 
from the NED database all the redshifts in a 1 deg 2 region 
centered on the cluster position. The histogram of these red- 
shifts shows two concentrations around z— 0.03 and z— 0.15. 
Since these concentrations are present over the 1 deg 2 field 
of view, we suggest that there are two foreground galaxy 
filaments or sheets. In this regard, Fukugita et al. (1995) 
predict i' — z' — 0.3 for early galaxy types at these (both 
— 0.03 and — 0.15) redshifts, which implies these filaments 
or sheets are populated by relatively evolved systems. 



2.2.3. Near IR observations 

We also obtained near infrared data (J and Ks bands) in 
order to fill the gap between the i' and z' Gemini bands and 
the Spitzer infrared data. 

A J band Subaru MOIRCS image was taken with a 
40 min exposure in 2007. Data were reduced with the stan- 
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Fig. 3. The color magnitude relation computed from i' and 
z'. The galaxies with spectroscopically determined redshifts 
that place them in the cluster (0.850 < z < 0.874) are 
red diamonds. Galaxies with spectroscopically determined 
redshifts outside the cluster are blue triangles. The black 
dots are all galaxies. The green squares mark the only two 
galaxies at the cluster redshift that were not detected in 
the MIPS 24 fim band. 




Fig. 4. Histograms showing the completeness levels for the 
i', z', J, and Ks bands. 



dard MOIRCS pipeline and provided a catalog complete 
down to ,1-23 (sec Fig.^). 

We also observed during 2 hours in the Ks band a sub- 
field of the RX J1257 area at the 3.5m APO telescope with 
the GRIM camera in 2000. This provided a catalog com- 
plete up to Ks— 22 (see Fig. |4j in a smaller area (see Fig. [5j). 



2.2.4. Catalog extraction 

The angular resolution of our images ranged from a few 
tenths of an arcsec (for the i' GMOS — 5/2 x 5/2 image) 
to a few arcsec (for the IRAC and MIPS images). It is 
therefore impossible simply to extract object catalogs in- 
dependently from all these images and to correlate these 
positions because of source confusion. We therefore chose 
to align and rebin to the same p ixel size (0.15 arcsec) a ll our 
images using the Terapix tools |http: / /terapix.iap.fr/[ ) and 
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Fig. 5. CMOS ~ 5/2 x 5/2 field i'-band plus the positions 
of the galaxies located in the z =[0.82,0.90] interval (de- 
rived from LePhare) over plotted . Black isocontours are 
computed at the positions of these galaxies with an adap- 
tive kernel technique and a smoothing window of 7 arcsec 2 . 
The fit templates have types from the LePhare web page 
(see footnote in § l2.1.3p . between and 30, being the ear- 
liest type and 30 the latest (star-burst) type. Red, types 
between and 10, green between 11 and 20, and blue be- 
tween 21 and 30. 



the standard procedure developed for the CFHTLS fields 
(McCracken et al. 2008). We used the weight maps for each 
band to produce two areas: one being the common area 
of all images and the other being the common area of all 
images except the Ks image (which is significantly smaller 
than all the other images). This is illustrated in Fig. [51 

We then extracted an object catalog from these images 
with the SExtractor package, in double-image mode, with a 
detection threshold of 1.1 and a minimum number of pixels 
(above this threshold) of 3. This provided an i', z', J, Ks, 
3.6/im, 4.5/im, 5.8/im, 8.0/im, and 24/^m catalog. 

2.2.5. Photometric redshift computations 

We have computed photometric redshifts with the LePhare 
package. A full description of this tool is given in Ilbert 
et al. (2006) and in the LePhare web page, but here are 
the salient points. LePhare is based on a template-fitting 
procedure. We used templates from Polletta et al. (2007) 
and star- forming templates from Bruzual & Chariot (2003). 
The zero-point calibration has been done on the basis of 
the spectroscopic catalog with an iterative procedure de- 
scribed in section 4.1 of Ilbert et al. (2006), by comparing 
the predicted magnitudes from the best fit template and 
the observed magnitudes. 

LePhare required shifts of 0.07 in i', —0.24 in z' (data 
observed under non-photometric conditions), 0.22 in J, 
—0.44 in Ks (data observed under non photometric con- 
ditions), 0.135 at 3.6 ^m, 0.08 at 4.5 /jm, 0.28 at 5.8 /j,m, 
0.22 at 8.0 //m, and 0.0 at 24 /im. 

The relatively high values of these shifts can be ex- 
plained in part by the correlation between the shift am- 
plitude and the photometric quality of the night. 

The shift in the photometry can also be explained in 
part by the nature of the spectral catalog that we used 
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Fig. 6. Spectroscopic versus photometric redshifts. Open 
symbols: multi-peak PDFs, where the redshift chosen is 
the one with highest probability; filled symbols: single-peak 
PDFs. See text. 



which is dominated by galaxies in the RX J1257 cluster. 
However, the spectral shapes of cluster galaxies are not well 
characterized in the only available LePhare template library 
which contains only field galaxy templates. Experience has 
shown that comparing field templates to cluster galaxies 
usually leads to high shift values. For example, similar shifts 
of 0.15 were applied to the Coma cluster photometric data 
by Adami et al. (2008). 

A set of fit parameters was then produced and we are 
mainly interested in multicolor type (based on a classifica- 
tion in a 9 magnitude-color space), redshift, and PDF es- 
timates. We first assess the quality of the photometrically 
determined z (z p , photo-z) by comparison with the spectro- 
scopically determined z (z s ), assuming that the uncertain- 
ties in z s are negligible. If Az — z s — z p , we can estimate 
the redshift accuracy from <jaz/(i+z) using the normalized 
median absolute deviation (NMAD Hoaglin et al. 1983) de- 
fined as 1.48x median(|z p — z s \/(l + z s )). This dispersion 
estimate is robust with respect to catastrophic errors (i.e. 
objects with \z p — z s \/(l + z s ) > 0.15). The percentage of 
catastrophic errors is denoted by r\. 

We show in Fig. [5] the relation between the estimated 
photometric redshift and the spectroscopic value of the red- 
shift. Beyond the general agreement, we see that the pho- 
tometric redshift of several galaxies is not very well deter- 
mined. Most of these galaxies have several peaks in their 
PDFs (open symbols in Fig. [6]) and this shows the need to 
consider the PDF instead of the best fit value of the pho- 
tometric redshift. I me dian is the median magnitude of the 
spectroscopic sample taken as reference. We used l/(l+z) 
weighting of the residuals to be consistent with previous 
work such as Ilbert et al. (2006). 

Based on the PDF, we computed the integrated proba- 
bility for each galaxy to have a redshift inside a given in- 
terval. From the spectroscopic catalog, we determined that 
the optimal interval is z=[0.82,0.90]. When predicting that 
a galaxy is within this interval on the basis of the PDF, our 
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success rate is close to 80%. When predicting that a galaxy 
is outside this interval, our success rate is close to 65%. 

We therefore extracted a candidate list of galaxies with 
a high probability of being part of the cluster by being 
included in the interval z = [0.82,0.90]. We plot in Fig.[5]the 
position of these galaxies along with their template color 
type based on the fits of the SEDs with LePhare. We note 
that the X-ray centroid is within the contours (see also 
Fig- HHdiscussed later in the text) of this concentration of 
galaxies. 

3. Discussion 

3.1. Identification as a cluster 

The first evidence that RX J 1257 is a massive cluster, is 
the strong peak at z ~ 0.87 seen in FigJT] 

Second, although there is a relatively uniform spatial 
distribution of galaxies at the cluster rcdshift, there is also 
a significant concentration at the putative cluster location 
as defined by the galaxy concentration contours (see Fig. [5]) 
and the X-ray emission (see Figs. [3] and [5]) ■ Thus, this sys- 
tem is likely a cluster embedded in a larger filament or 
sheet of galaxies, similar to the sheet of galaxies detected 
at - 0.735 in the CDFS (e.g. Adami et al. 2005). 

Third, we applied the SG (Serna & Gerbal 1996) 
method to our rcdshift catalog. This hierarchical method 
allows substructures to be extracted from a catalog con- 
taining positions, magnitudes, and rcdshifts, based on the 
calculation of their binding energies. It provides a robust ve- 
locity dispersion and a virial mass. The SG method reveals 
the existence of a group of 18 galaxies located at the can- 
didate cluster center, coincident with the X-ray emission, 
with a mass of 6.1 x 10 14 M Q , based on a velocity dispersion 
of 600 km s . Within this group, we found two subgroups 
of 5 and 3 galaxies with respective masses of ~ 10 14 M Q 
(velocity dispersion of 289 km s^ 1 ) and -2.2 10 13 M Q (ve- 
locity dispersion of 255 km s _1 ). These masses from galaxy 
velocity dispersion measurements are consistent with the 
inferred X-ray mass of (1 — 5) x 10 14 M Q , and are in the 
range to call RX J1257 a massive cluster (e.g. Sarazin 1986). 

Fourth, the X-ray emission is clearly extended, as seen 
in Figs. [1] and [2] and is similar in size to the galaxy concen- 
tration shown in Fig. [5J 

Fifth, assuming a rcdshift of 0.866, the projected size of 
the X-ray emission is about 0.5 Mpc, and its kT (3.6 keV) 
and inferred Lx.boi (2 x 10 44 ergs s _1 ) are all consistent with 
a massive cluster at 0.866 (see, e.g., Ettori et al. 2004a, and 

Fig. USD- 

3.2. Is RX J1257 bimodal? 

The XMM-Newton data suggest the X-ray emission is bi- 
modal, but the Chandra data show there is a point source 
located in the western portion of the X-ray emission. We 
show in Fig. [5] a smoothed (by an 11 arcsec radius Gaussian 
kernel) low energy (0.3-2.0 keV) Chandra image with XMM 
contours (red) superposed. There are clearly two centroids 
enclosed by the XMM-Newton contours (red in Fig. [H]) and 
the west-most centroid has a round (point) source near the 
center. This point source is probably the major contribu- 
tion to the western peak in the XMM-Newton contours as- 
sociated with the cluster. Nevertheless, as can be seen in 
Fig. [5] it appears that diffuse emission surrounds this point 
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Fig. 7. Histogram of the GMOS redshifts, with the region 
around 0.87 expanded in the insert. 



Fig. 8. A 11 arcsecond Gaussian smoothed Chandra image 
in the 0.3-2.0 keV range. The cyan contours are based on 
this Chandra image and arc drawn to highlight the Chandra 
data. The red contours are based on the XMM-Newton 
0.3—6 keV image (see §3.2.) and they coincide with those 
in green or yellow in Fig. [5J The text in the figure points 
to the two brightest regions in the Chandra image that are 
centered in XMM-Newton contours. 



source. If the point source were absent, the surrounding dif- 
fuse emission could be the origin of the second peak in the 
XMM-Newton contours. We conclude therefore, that there 
is still a centroid in the western X-ray contours that is pos- 
sibly not due to the point source (the nature of which is 
discussed below), and that the X-ray emission associated 
with RX J1257 has bimodal characteristics. 

This hypothesis of bi-modality fits well with our finding 
that the distribution of galaxy cluster members is also bi- 
modal (based on the SG method) with the two subgroups 
also aligned in the E-W direction. That the western sub- 
group of galaxies lies outside the XMM-Newton contours 
may indicate the relatively low dynamical age of the cluster 
(see Fig. [S]). This is confirmed by Fig. [5] where we see that 
late type galaxies are relatively numerous in RX J1257, be- 
ing about half of the cluster population in the central area. 

3.2.1. Nature of the embedded point source 

The point source shown in Fig. [TU] and located at 
a=194.2930, (5=47.6339 (from the Chandra data) could be 
associated with an optical object located at the approx- 
imate coordinates a=194.2937, 5=47.6340. However, due 
to the proximity in projection of three objects, SExtractor 
provided only one detection. With this caveat in mind, we 
estimate that the probability of this object to be a cluster 
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Fig. 9. A color image based on Spitzer IRAC 3.6 /jm (R), Gemini z' (G) and Gemini i' ( B) image. The XMM-Newton 
contours in the 0.3-6 keV interval are shown in green, left, and yellow/ white, right (adjusted to separate the northern 
AGN from the cluster emission). The green circles show galaxies at the cluster redshift. The magenta and blue circles 
indicate galaxies in the two sub-groups detected by the SG (Serna & Gerbal 1996) method; only 3 of the 5 objects in 
the western sub-group are shown here. Big and small white circles indicate the presence or absence of [Oil] emission: no 
white circles = no detectable [Oil], small white circle (only two, one in the southeast corner of the image, the other just 
north of the guide star mask) = [Oil] smaller than 2 times the noise, large white circles = [Oil] stronger than 2 times 
the noise. The rightmost image is a blow up of the left one. North is up and East to the left; a concordant cosmology 
is assumed, thus 463 kpc = 1 arcmin (as depicted by the horizontal white bars in both images) at the redshift of RX 
J125712+4738.17 (z = 0.866). 



Fig. 10. i'-band image of three blended galaxies that run 
from North to South in the image. The two lower most 
probable positions are marked by white circles. The lowest 
declination object was detected in the Chandra image with 
probably diffuse emission around it (see text). The next 
object to the North is most likely a MIPS source which is 
the 7th brightest MIPS object in the field. This region can 
be found in the northernmost of the two magenta circles 
that overlap in Fig [9j North is up and East is to the left 



member is 89tQ- The nature of this object (normal galaxy or 
active object) remains unclear, however, because a galaxy 
template provides a better fit to the composite spectrum 
than a QSO template, but the X-ray flux is 4.8 x 10~ 13 
ergs cm -2 s^ 1 which at the cluster distance means its Lx 
is about 10 45 ergs s _1 , typical of QSOs. The most con- 
servative explanation, then, is that this X-ray source is a 
chance alignment with the cluster and the derived photo-z 
(and spectral classification) for the merged (by SExtractor) 
source does not apply to the X-ray source. 
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Fig. 11. Histogram of the distribution of the magnitudes 
of the galaxies for which spectra were taken. The cross 
hatched area is for those galaxies for which redshifts were 
successfully derived. 



4 This was done by integrating the redshift probability distri- 
bution given by LePhare over the total redshift range. We then 
integrated over an estimated redshift range for the cluster of 
0.8-0.9, and calculated the ratio of these two values. See Adami 
et al. (2008). 



3.3. Characteristics of the cluster galaxies 

3.3.1. Selection Effects 

In Fig.[TT]we show the brightness distribution of the galax- 
ies for which we attempted to determine redshifts. ^From 
this, it can be seen that we were successful a very large 
fraction of the time and that there was no initial choice 
that would have favored emission line galaxies. Also, it was 
possible to determine redshifts even when emission lines 
were not present (over 50% of the time). Thus, we conclude 
that the relatively large faction of emission line galaxies we 
found at the cluster redshift was not due to a selection ef- 
fect, and we proceed in the following discussion based on 
this assumption. 

3.3.2. Spectra 

We first consider those galaxies with spectroscopically de- 
termined redshifts within the putative cluster redshift range 
(0.850-0.874). We find that about 45% (10/21) of the cluster 
galaxies are emission line objects (cf. Fig. [§]). These results, 
although based on small number statistics, are consistent 
with the hypothesis that RX J 1257 is relatively young and 
there has not been time for the cluster environment to cause 
the galaxies to change from late type (star forming) to early 
type (passively evolving). In comparison, in more mature 
low z clusters, ~75% of galaxies can be classified as early 
type objects (e.g. Adami et al. 1998). 

The fraction of MIPs detections within the population 
of spectroscopically confirmed galaxy members is nearly 
100% (19/21). This high fraction of MIPs detections could 
be explained if galaxy harassment has triggered star for- 
mation or if the system is so young that dust has not yet 
been swept from the late type galaxies. That the large ma- 
jority of the MIPS cluster members in the core are associ- 
ated with the western (rather than eastern) portion of the 
cluster X-ray emission (see Fig. [T2~|) is suggestive that the 
western half is the younger part of the cluster. The data 
also suggest that the system to the west is in the process of 
merging with the group to the east. Therefore, an appeal- 
ing scenario that explains our results is that proposed by 
Elbaz et al. (2007) where late type galaxies at z~0.9 are red 
due to dust. These galaxies do not become blue until they 
have lost enough dust, which could come about through a 
variety of evolutionary processes. 

3.4. Cluster Collapse 

We suggest this cluster is just in the process of formation 
based on: (a) the bimodal distribution of both the X-ray 
emission and the galaxy population; (b) the fact that the 
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Fig. 12. Left: a MIPS 24yum (R) , IRAC-3.6 (G) and i'-band (B) image that shows all the spectroscopically determined 
members that were detected by the MIPS (red boxes; based on MlPS-i-band source association that is not necessarily 
correct to better than 2 arc seconds) with our observations. The circle is 3 arcmin radius or about 1.3 R v (see text). 
Right: the same color scheme on a zoomed scale with the i-band (blue) now greatly suppressed to bring out the MIPs 
flux. The green contours are based on the 0.3-2 keV XMM-Newton image. North is up and East is left 



majority of the spectroscopically confirmed cluster mem- 
bers were detectable in the MIPS 24 /im channel; (c) the ex- 
istence of cluster galaxy substructures based on the Serna- 
Gerbal method; (d) the kT on the high side relative to the 
predicted Lx,boi-kT relation (the thin dotted line in Fig. I13|) . 

If there is some energy input prior to infall, this would 
produce a higher kT versus Lx,boi than we would expect 
from energy input from the infall alone. For example, the 
initial energy injection in addition to gravitational infall 
comes from some other process such as galaxy outflows or 
SN (see Babul et al. 2002, Jones ct al. 2003, Bode et al. 
2007). Since we have found several late- type cluster mem- 
ber galaxies, this suggests that galactic outflows have not 
dominated the heating of the RX J1257 ICM, otherwise gas 
and dust loss should have made these galaxies faint rather 
than bright 24 /im objects. 



lations that showed that clumpy, not fully merged clusters 
will tend to have luminosities that fall below the Lx,boi-kT 
relation for a given kT. Thus, the X-ray data are consistent 
with our expectations that the cluster is in the process of 
merging or collapsing and is not fully relaxed. 

3.5. A Comprehensive Picture 

The results presented here combined with previous work 
(references above) on clusters of galaxies lead to a scenario 
in which part of the red galaxies in high z clusters are red 
due to dust rather than being red and dead. The less mas- 
sive the cluster, the younger it will be in terms of having not 
yet completed infall, compared to more massive clusters. A 
system such as RX J1257 should then, as we found here, 
have a high fraction of red dusty galaxies. These galaxies 
will become blue after dust ejection and then red again after 
having consumed all their gas. 



100.0 



10.0 



1.0 



0.1 



+ 0.4 < z < 0.6 
x 0.6 < z < 0.8 

x 0.8 < z < 1.0 
A z > 1.0 




10 



(kcV) 



Fig. 13. Based on the Ettori et al. (2004a) Fig. 4, right 
upper panel, with the data point added for RX J1257 as a 
green cross. This is the L — T relation with correction by 
E(z). Dotted line: slope fixed to the predicted value of 2. 
Dashed line: slope free. See Ettori et al. (2004a) for further 
details. 



Another way to consider Fig. [T3] is to assume a lower 
than normal Lx,boi for a given kT. Such a low Lx,boi is ex- 
pected by Ventimiglia et al. (2008), who carried out calcu- 



4. Summary and Conclusions 

We have discovered the young massive cluster of galaxies 
RX J1257.2+4738, at z = 0.866. The cluster has a lumi- 
nosity of about 2 x 10 44 ergs s _1 , a kT of about 3.6 keV, 
and an estimated mass of 6 x 10 14 M Q . 

The cluster RX J1257 is in the process of collapsing. At 
birth, there was possibly significant energy input to the in- 
tracluster medium besides infall which produced a cluster 
having an intracluster medium temperature which is not 
the one expected for its luminosity. The MIPS 24/im detec- 
tions plus [Oil] 3727 emission line galaxies imply a high star 
formation rate in several of the cluster member galaxies. 

In agreement with Elbaz ct al. (2007), we suggest that 
there should be an evolutionary sequence from red dusty 
galaxies in a cluster to blue galaxies and finally to red-dead 
galaxies. 
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